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SUMMARY OF FLUTTER EXPERIENCES AS A GUIDE TO 
THE PRELIMINARY DESIGN OF LIFTING 
SURFACES ON MISSILES 1 
By Dermis J. Martin 

SUMMARY 


A limited review is made of some experiences in the flight testing 
of missiles and of wing flutter investigations that may be of interest 
in missile design. Several types of flutter which may be of concern in 
missile studies are briefly described. Crude criteria are presented 
for two of the most common types of flutter to permit a rapid estimate 
to be made of the probability of the occurrence of flutter. Many of the 
details of the flutter picture have been omitted, and only the broader 
elements have been retained so as to give the designer an overall view 
of the subject. 


INTRODUCTION 


Many different types of flutter may be encountered on airplanes , 
propellers, helicopters, and missiles and the speed ranges and conditions 
encountered lead to flutter phenomena that are widely different. Broadly 
speaking, the phenomenon of flutter is generally concerned with vibra- 
tions or oscillations of a lifting surface. Oscillations of a lifting 
surface give rise to oscillations of the aerodynamic forces which in 
turn, under certain conditions, may have phase characteristics that 
increase the oscillations to dangerous amplitudes. Some types of flutter 
may be mild; others may be disastrous. Flutter may involve fully estab- 
lished flow or broken-down flow, high or low frequencies of the structure, 
and one or more modes of vibration. 

The missile not only experiences many of the flutter problems 
encountered with airplanes but also presents many new and different 
problems, depending upon the design and purpose of the missile. Examples 
are: skin flutter, flutter of automatic controls or servomechanisms, 

-^-Supersedes declassified NACA Research Memorandum L5U50 by 
Dennis J. Martin, 1951- • 
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and flutter of short wings with ram jets or external stores. Many of 
these types of flutter can best be studied by difficult experiments; 
others require long and tedious theortical investigations. For the more 
common types there exist sufficient experimental data to evaluate simple 
criteria. In general so many factors enter into a flutter case that a 
comprehensive criterion becomes quite unwieldy. Simple criteria must 
neglect or restrict many parameters. Furthermore, there are possibilities 
for exception; hence any simple criterion should not be considered as 
perfectly general. In spite of these limitations a criterion does have 
some usefulness in estimating the probability of a particular type of 
flutter occurring for a given configuration. 

In this paper two simple criteria are presented. The first is for 
the most common type, the wing bending-torsion flutter. Another is 
presented for stall flutter, and a brief discussion is included of pitch- 
bending flutter. 


SYMBOLS 


A 


panel aspect ratio 


a 


velocity of sound; also, nondimens ional distance of elastic 
axis behind the midchord expressed as a fraction of the 
semi chord 


b 


semichord 


c 


chord 


f^A) approximate function relating torsional frequency of a tapered 

panel to that of an untapered panel 


f 2 (A) 


function relating the chord of a tapered panel at the panel 
three-quarter span to the chord at the panel midspan 



effective shear modulus of an equivalent section 


I 


moment of inertia about elastic axis 


J 


section torsional modulus 


J, 


S 


section torsional modulus of a solid section 


JG 


torsional stiffness parameter 


* * 
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i L 

panel span 

M 

Mach number 

m 

mass per unit span 

N 

calculated parameter (two-dimensional incompressible flutter 
speed divided by velocity of sound) 

n 

empirical number used in aspect-ratio correction 

P 

fluid pressure 

r a 

nondimens ional radius of gyration relative to elastic axis, 
expressed as fraction of semichord 

t 

section thickness 

V 

velocity 

v f 

& 

flutter velocity 

X 

nondimens ional geometric parameter, defined in appendix 

X 

distance of wing leading edge behind the missile center of 
gravity 

x a 

nondimens ional distance of section center of gravity behind 
elastic axis expressed as fraction of semichord 

a 

angle of attack measured from zero lift 

A 

sweep angle of wing midchord line 


taper ratio (ratio of tip chord to root chord) 

€ 

nondimens ional distance of section center of gravity behind 
wing quarter -chord expressed as fraction of chord 

K 

relative density parameter 

P 

fluid density 


angular wing torsional frequency 


angular wing bending frequency 


Subscripts : 


0.50 


panel midspan station 
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0.75 panel three-quarter-span station 

0 sea-level standard conditions 


DISCUSSION OF CRITERION 


Types of Flutter 

Examples of some of the more common modes that may interact during 
flutter are given in figure 1. The first example shown is the most 
common type of flutter encountered; in this type of flutter the elastic 
modes of the wing (wing bending and wing twisting or torsion) combine 
to extract energy from the air stream, that is, to produce flutter. In 
addition, a control surface may interact significantly with these motions 
to produce other types of flutter. The second example shown in figure 1 
illustrates a type of flutter which involves only one motion or degree 
of freedom. The type of flutter illustrated occurs at high angles of 
attack and is commonly known as stall flutter. Only a torsional twisting 
motion of the wing is present. There are other motions that may produce 
a single -degree flutter of the type illustrated by this stall-flutter 
case. Examples are: aileron buzz, single-degree bending oscillations 

of swept wings, and single-degree pitching oscillations of a wing. The 
third example in figure 1 illustrates a type of flutter in which the 
motion of the whole fuselage enters significantly into the flutter. 

This example illustrates a pitching motion of the entire missile combined 
with a bending motion of the wing. Other body motions (rolling, yawing, 
vertical translation, and so forth) also may enter into flutter. 


Bending -Torsion Flutter 

In order to illustrate the significance of the first criterion to 
be presented for the most common flutter, wing bending-torsion flutter, 
figure 2 (see references 1 and 2) has been prepared to show the flutter 
behavior of wings over a range of Mach numbers . Shown in this figure 
is the actual flutter Mach number plotted against a calculated quantity 
N, which is dependent upon the wing stiffness, center-of -gravity loca- 
tion, mass ratio, certain aerodynamic quantities, and so forth. The 
parameter N is an orderly combination of many of the parameters that 
are important in flutter. The flutter expert may recognize this N 
as the calculated two-dimensional incompressible flutter speed divided 
by the velocity of sound. Topical curves for wings of zero sweep and 
full-span aspect ratios of 2 and 7 are shown, and an approximate curve 
is shown for a wing of 60° sweepback and aspect ratio 4. As the value 
of N for the aspect-ratio-7 wings is increased, for example, by 
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increasing the wing torsional stiffness or by increasing the operating 
altitude, the flutter Mach number is seen to increase and a value of N 
is eventually reached which will not produce an intersection with the 
flutter boundary. The significance of this result is that, if the wing 
is so stiff or if the quantity N is so large that flutter is not encoun- 
tered in this region around a Mach number of 1 (for the higher aspect 
ratios), flutter of the type considered herein would not be expected 
to occur at higher Mach numbers . It is this maximum or critical value 
of the quantity N that is of interest in missile design, because mis- 
siles must operate throughout the Mach number range. It must be empha- 
sized that the curves for other aspect ratios, sweep angles, and so forth 
may appear quite different from the one shown for the aspect -ratio -7 wings 
and hence may have a different critical Mach number range and a different 
value of N required for the wing to be flutter-free. If all the criti- 
cal values of this parameter N that are necessary to avoid flutter were 
known for the various aspect ratios, sweep angles, thickness ratios, 
section properties, and so forth, the flutter problem for missile design 
would be greatly simplified. 

Investigations in the Langley 9- by 18-inch supersonic flutter tun- 
nel have attempted to define values of N for various supersonic Mach 
numbers. (See refs. 3, 4, and 5 .) Many of these data, together with 
many data from the Langley Pilotless Aircraft Research Division on mis- 
siles and bombs which experienced no known flutter difficulties, as well 
as on missiles on which flutter was attained, have been accumulated. 

This store of experience has been compared to a simplified and modified 
criterion which groups the significant parameters in a manner similar to 
that used for N in an attempt to establish limits of the critical values 
of the structural and aerodynamic requirements for a wing to be flutter - 
free. This criterion is based on modifications to an approximate flutter 
formula proposed by Theodorsen and Garrick (ref. 6 ) . This formula was 
for high-aspect -ratio, heavy wings having a low ratio of bending to tor- 
sional frequency. The application of modifications of this formula to 
include low-aspect -ratio wings including swept and highly tapered wings 
is admittedly stretching the basic formula; nevertheless, in the present 
study this approach has been made, and the parameters have been adjusted 
■until a reasonable coherence in the results was obtained. For simplicity 
this modified criterion is broken down into simple geometric dimensions 
and structural properties. (This modification is described in the 
appendix.) The experimental data which have been accumulated are then 
compared, and an attempt is made to bracket the safe wings and the unsafe 
wings. This comparison is made in figure 3 . 

Plotted against the effective shear modulus of the wing material 



2 


6 


NACA TN 4197 


obtained from the geometric dimensions of one wing panel. (See appendix.) 
The abscissa G™ is the effective shear modulus of the wing structure 

Cj 

and is indicated for wings of solid wood, magnesium, aluminum, titanium, 
and steel. A solid wing of, say, aluminum would fall at the point marked 
along the abscissa while a fabricated wing of aluminum would have a 
lower effective G E and would fall somewhat to the left, depending upon 

the skin thickness and spar size. The value of G^ can be determined 

for fabricated wings from a measured value of the torsional stiffness 

parameter JG by the relation Gg = Treasured, -where c is the 

ct^ 

chord and t is the thickness. The quantity X noted in the ordinate 
is shown in figure 4 as a function of the panel aspect ratio for constant 
values of the thickness ratio in the streamwise direction. It must be 
remembered that the abscissa of this figure is the exposed aspect ratio 
of only one wing panel as distinguished from the normal aspect ratio 
which includes both wings and the fuselage . 


Shown in figure 3 are data taken from subsonic, transonic, and 
supersonic wind tunnels and from rocket and bomb-drop tests for both 
swept and unswept wings. The open points are for missile wings that 
traversed the Mach number range to at least a Mach number of 1.3 or 
higher without known failure. The solid points are for missile and. 
wind-tunnel tests where flutter or failure occurred. It must be pointed 
out that some of the data are for missiles that were designed primarily 
for aerodynamic research. The instrumentation of these missiles was, 
therefore, not usually of a type that could definitely indicate that no 
oscillations occurred for the cases represented by open points or that 
the failures for the cases represented by the solid points were due to 
flutter rather than some other cause. The many data shown tend to indi- 
cate that two regions can be defined in which the open and solid points 
-are reasonably well separated and the flutter region is established. 

The shaded area indicates a probable division based on the existing 
data. This chart is useful in estimating the probability of the occur- 
rence of flutter of the bending-torsion type for a given configuration. 

A designer may see where a given design lies with respect to many other 
designs which did or did not experience flutter problems . 


As an illustration, if a design had an exposed-wing-panel aspect 

ratio of 2 and a streamwise thickness ratio of 4 percent, from fig- 

c 

ure 4 a value of X of about 1.23 x 10 pounds per square inch is 
indicated. If the missile were ground-launched, that is, at standard 

pressure -E— = 1, and if the wing were untapered, that is — — 1, 

^0 5 

then the ordinate of figure 3 for this design would be 1.25 x 10 ; and 

if the construction were solid magnesium, it would plot in the flutter 
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region and would most probably be unsafe. If the wing were of solid 
aluminum, it would be marginal. Further detailed analysis or experiment 
would be needed to complete this design. However, if the wing were of 
solid steel, it would probably be safe, at least insofar as the bending- 
torsion type of flutter considered is concerned. 

As another example, suppose the material of construction has already 
been selected, titanium for instance. In order to allow a reasonable 
margin of safety an ordinate of figure 3 of about 0.8 or less might be 
specified. If the wing were untapered and ground - launc hed , the ordinate 
of figure 4 is then 0.8 and it can be seen that for a panel aspect ratio 
of 1 a thickness ratio of 2-5 percent is required. For an aspect ratio 
of 2 a thickness of 4.5 percent is required, and for an aspect ratio of 3 
a thickness of about 6.5 percent must be used so that the design may 
most likely be free of the bending-torsion type of flutter for which 
this f igure applies . 


Stall Flutter 

There are many other types of flutter that may occur under certain 
conditions, and they must also be investigated. With a change in the 
type of flutter, a change must be made in the type of criterion. As 
mentioned previously the type of flutter that may be encountered depends 
upon the design and purpose of the missile. As an example of the depend- 
ence of the type of flutter upon the use of the missile, it may be men- 
tioned that high-angle -of -attack flutter, that is, stall flutter, would 
probably be considered as possible only for missiles that are required 
to maneuver sharply . This conclusion is for ground -launched missiles } 
however, any air-launched missile that is carried externally may be sub- 
ject to large angles of attack during airplane maneuvers prior to 
launching and thus may become subject to stall flutter. During a recent 
bomb-drop test at Langley, a missile-wing failure occurred while the bomb 
was attached to the airplane. The failure occurred at speeds consider- 
ably below the flutter speed subsequently obtained with an identical wing 
that was protected from the airstream while attached beneath the 
airplane . 

An investigation of stall flutter of thin wings was therefore begun. 
Although data on stall flutter encountered on missiles are not readily 
available, a brief discussion of the stall flutter of thin wings and 
stall flutter of propellers may serve as a rough guide for missile design. 

Figure 5 illustrates the flutter behavior of a typical wing at low 
speeds as the wing produces lift. The ordinate is a nondimens ional 
flutter-speed coefficient V/ bo> a , where V is the flutter speed, b is 

the half -chord, and o> a is the torsional frequency. The abscissa is 
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the angle of attack. As the angle of attack is increased, the flutter 
speed is reduced drastically. The flutter speed falls rapidly and a 
minimum is reached at an angle of attack near the stalling angle of the 
wing. 


The flutter encountered at the low angles of attack is of the 
bending-torsion type which was discussed in the previous figures and 
was shown to be strongly dependent upon the material of construction. 

At high angles of attack the flutter occurs essentially in only a tor- 
sional mode, and this minimum value of the flutter-speed coefficient 
V/bo^ has been found to be nearly equal to 1 for almost all wings and 
propeller blades at low speeds, regardless of the material of construc- 
tion. These results have been confirmed for both wings and propeller 
blades, and the results are thought to be generally valid for the 
subcritical-f low speed range . What these curves might look like at super- 
sonic speeds has not been determined; however, a study of the minimum 
value of V/bo^ as affected by Mach number has indicated a beneficial 

effect at higher speeds . The experimental work of Baker at Langley 
(reference 7) has suggested that the quantity boo a referred to the 

speed of sound was a significant parameter for determination of con- 
figurations that would be free of stall flutter . 

Figure 6 has been prepared to show a comparison of experiment with 
this parameter bcjc^/a for a range of Mach numbers. Shown in this fig- 
ure is the flutter Mach number plotted against bcjo a /a. The curve shown 

represents the boundary where stall flutter could begin for a given 
value of bco a /a. As bcjo a /a is increased, for example, by increasing 

the chord or increasing the torsional frequency, a value is noted to be 
reached which will not produce an intersection with the flutter boundary. 
The result is quite similar to the situation that occurred for the 
bending-torsion flutter (fig. 2). 

Baker (ref. 7) has shown that a value of boe^/a of at least 0.5 

is required for a propeller to be completely free of stall flutter. 

Rainey (ref. 8) has substantiated this value of 0.5 for unswept wings 
of moderate aspect ratio and low structural damping but has indicated, 
however, that aspect ratio, structural damping, and sweepback may have 
an influence on the critical value of b<u a /a. These effects are not 

well-determined and cannot at present be included in a design chart. 

The value of bcq^/a of 0-5 has nevertheless been used to prepare a 
design chart for solid unswept wings . Since the torsional frequency 
times the chord for solid unswept wings is a function only of the length- 
chord ratio L/c and the thickness ratio and is essentially independent 
of the material (that is, for such common materials of construction such 
as steel, aluminum, and magnesium), the design chart (fig. 7) is presented 
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in terms of the length-chord ratio and the thickness ratio required to 
attain a value of bco a ^a of 0-5* A speed of sound of 1,100 feet per 

second is assumed. 

Wings having geometric quantities which plot above the solid line 
of figure 7 may experience stall flutter if sufficiently high angles of 
attack are encountered. It must be remembered that the boundary indi- 
cated is for the most critical condition. If the speed range is trav- 
ersed at low angles of attack, a design may be well above the boundary 
of this figure without encountering difficulties. This boundary repre- 
sents only the conditions required for the wing to be completely free 
of stall flutter throughout the speed range at any angle of attack. 

The margins of safety for this criterion are not established, and the 
criterion may have to be modified as more information and data are made 
available on sweepback, aspect ratio, and structural damping. 


The significance of free-body modes in flutter has been of interest 
for some time. The problem was considered in early British work on 
flutter involving the mobility of the fuselage. For example, Broadbent 
(reference 9) developed simple criteria based on the position of the 
nodal line. The type of flutter involving missile pitching and wing 
bending is dependent upon the moment of inertia in pitch of the missile 
and the bending stiffness of the wing as well as upon the wing location 
with respect to the center of gravity of the missile. For wings which 
meet the torsional-stiffness criterion for the bending -torsion type of 
flutter but are weak in bending, this type of flutter may become impor- 
tant for some wing locations if the missile has a high moment of inertia 
in pitch. In some flutter tests of rocket vehicles at Langley, several 
failures have occurred and seemed to involve principally wing bending 
and missile pitching. The frequency of flutter was somewhat below the 
first-bending frequency of the wing, near the short-period oscillation 
of the body. Analyses have been made (ref. 10), and the effect of wing 
location is illustrated in figure 8. 


in this case the first-bending frequency is used. The abscissa 

is the nondimens ional distance of the wing behind the center of gravity 
of the body. The dashed line represents the conventional bending- 
torsion type of flutter while the solid line shows the effect of inclu- 
sion of a body degree of freedom. This flutter speed is much lower than 
the bending-torsion type for rearward locations and much higher for for- 
ward locations. The significant conclusion that can be drawn from these 
studies of pitch-bending flutter is that the most important considera- 
tion is the inclusion of the proper degrees of freedom or modes in the 


Pitch-Bending Flutter 


The ordinate in figure 8 is the flutter-speed coefficient 
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analysis. Moreover, the observation can be made that, with the change 
in the type of flutter, a change occurs in the type of flutter criterion 
thus in this case the critical speed is affected strongly by the bending 
stiffness and not by the torsional stiffness as in the case of bending- 
torsion flutter. 


CONCLUDING REMARKS 


Simple flutter criteria have been presented to serve as a guide in 
the preliminary design of lifting surfaces on missiles. The proximity 
of a configuration to a bending-torsion type of flutter instability may 
be indicated, and estimates may be made of the probability of the occur- 
rence of flutter by comparison of a given configuration with other 
designs that have or have not experienced flutter. A simple criterion 
is presented for stall flutter, and although the margins of safety can- 
not be established, the criterion may be useful in preliminary design. 
Pitch-bending flutter is briefly discussed and an example is cited which 
shows that the wing-bending stiffness and the configuration center-of- 
gravity position may strongly influence the flutter speed. The discus- 
sion of pitch-bending flutter illustrates the case in which a change in 
the type of flutter can bring about a change in the type of criterion 
that is needed. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October l6, 1957* 
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APPENDIX 

DERIVATION OF SIMPLIFIED FLUTTER CRITERION 


An empirical expression for flutter speed as given by Theodorsen 
and Garrick (ref. 6) for heavy, high-aspect -ratio wings with a low ratio 
of bending frequency to torsional frequency is 


V f _ 2V f 
Ml 



i/2 

\ + a + % 


( 1 ) 


i + a + *a , 

If e is substituted for (the nondimens ional distance of 

2 

the center of gravity of the section behind the quarter -chord position) , 
equation (l) , when squared, becomes 


Vf 2 


c^d 2 r ^ 
1 a 

16 k € 


( 2 ) 


Experience has generally indicated that more realistic values of 

are obtained if the quantities in equation ( 2 ) are evaluated at the 
three-quarter spanwise station. Geometrically similar sections are 
assumed at all spanwise stations; thus, equation ( 2 ) becomes 


v f 2 = c o.79^q 2r a 2 
a 2 l6itea 2 


(3) 


where a is the velocity of sound; also 



( 4 ) 
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and 


rtpb 2 _ flpc 2 

m 4m 


(5) 


The expression for the torsional frequency of a uniform beam is 



JG 

I 


(6) 


Coleman (ref. ll) has given the frequency of a tapered beam of constant 
thickness ratio in terms of the frequency of an untapered, uniform beam 
having the same root chord as 



tapered 


(a) f-,(A) 

\ a '/ untapered 


(7) 


where it has been found that f ^ 


may be approximated by 


f]_ w 1 + 1.87(1 - A ) 1-6 ( 8 ) 

When equations (4), (5), ( 6 ), and (7) are substituted into equation ( 3 ), 
the following is obtained: 



JtC Q.73 f 
4epc 1 < ~L 2 



(9) 


For a solid, thin airfoil J can be closely approximated by 



( 10 ) 
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For a fabricated section, J is extremely difficult to calculate. The 
value of JG of a section may be experimentally measured and an equiva- 
lent solid section may be assumed; thus 


(JG) 


measured 


j s g e 


( 11 ) 


where is an effective shear modulus of the equivalent section; 

that is, 

Ge = (12) 

ct' 5 


Since — is the thickness ratio and 
c 

equation (9) may be reduced by using 
relations : 


is the panel aspect-ratio, 

c 0.50 

equation (12) and the following 


^.75 ~ c o.50 f 2^^ 


(13) 


f 2 (A) 


i + y\ 

2(1 + A) 


(14) 


Equation (9) becomes, for wings 


of constant thickness ratio, 



I 2 = 5X1 

(i) 

'Vl 2f 2 2 

V a ) 

24epA^a^ 


(15) 


An aspect-ratio correction of the form has frequently been used. 

A + n 

Reference 12 has suggested a value of n = 2. This value of n = 2 has 
been used with success in flutter criteria in reference 13 . In the pres- 
ent investigation values of n = 1, 2, and 3 were tried and the value 


Ik 


NACA T N 4197 


of n = 2 gave the more consistent results. If the value of n = 2 
is used, equation ( 15 ) becomes 



( 16 ) 


It has been found that — — is closely approximated by ^ 1 

f 2 f 2 2 

1 2 

By use of the gas lav, the following is obtained: 


pa 


2 



(17) 


The ratio of specific heats 7 = 1*4 is used, and the results apply to 
air. If a different gas is used, the results must be corrected for a 
different specific-heat ratio. A value of € = 0.25 is assumed; how- 
ever, for sections with the center of gravity far from the 50 ~P ercen t- 
chord position a correction may be required. Equation (l6) becomes 



(18) 


Thus, if a critical value of exists, a plot of the denominator 

of equation (l8) for various materials or values of the numerator, Gg, 


for missile and wind-tunnel tests may permit a systematic separation to 
be made of the safe wings and the unsafe wings. The first portion of 
the denominator of equation (l8) 
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X = 


39-3A 3 

(i) ? (A + 2) 


(19) 


is calculated and plotted in figure 4. 

The value of X can be determined from figure 4 from the thickness 
ratio and aspect ratio of one exposed wing panel. This value of X is 

multiplied by — and by — to obtain the ordinate of figure 3. 

^ P 0 
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WING BENDING-TORSION FLUTTER 



Figure 1.- Examples of flutter modes. 




Figure 2.- Trend study of swept and unswept wings at transonic speeds 
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Figure 3*- Composite chart for bending- tors ion flutter. Ge 
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Figure 4.- The parameter X as a function of panel aspect ratio and 

thickness ratio. 
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Figure 5-- Flutter behavior of a typical wing at angles of attack. 



Figure 6.- Trend study of stall flutter. 
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Figure 7-- Composite chart for stall flutter of solid unswept wings. 


v/bo> h 



Figure 8.- Study of flutter involving body modes. 
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